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Abstract 

Exact quantum integrability is established for a class of multi-chain electron models 
with correlated hopping and spin models with interchain interactions, by constructing the 
related Lax operators and i?-matrices through twisting and gauge transformations. Exact 
solution of the eigenvalue problem for commuting conserved quantities of such systems is 
achieved through algebraic Bethe ansatz, on the examples of Hubbard and t — J models 
with correlated hopping. Our systematic construction identifies the integrable subclass of 
such known solvable models and also generates new systems including the generalized t — J 
models. At the same time it makes proper correction to a well known model and resolves 
recent controversies regarding the equivalence and solvability of some known models. 
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1 Introduction 

There is an upsurge of interest in recent years in the study of correlated electron and spin sys- 
tems in low dimensions, which is motivated mostly by the recent possibility of their fabrication 
and experimental verification of related theoretical results ]IJ, as well as by their potential for 
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applications to the high T c -superconductivity ||. On the other hand one can take an exclu- 
sive advantage in one-dimension, since some class of models in this case may become exactly 
solvable and even completely integrable, which therefore can enhance our understanding of the 
low-dimensional physics by providing detailed picture of the system through exact results. Few 
well known examples of such models are the XXX and XX Z spin chains [0], the Hubbard 
model [§], [|, the t — J model @, [7| etc., along with their various extensions and the ladder or 
their multi-chain generalisations P, j9[ fTOfl. Recently, a large class of models describing elec- 
trons with correlated hopping, spin ladders with interchain interactions and their multi-chain 
generalizations have attracted special attention due to their simplicity and exact solvability 
I!], [12], O, ITJL ITH]. A systematic analysis has also been carried out to identify the general 



class of these models which exhibits exact solvability by coordinate Bethe ansatz (CBA). How- 
ever, the important question of complete quantum integrability for such systems, which defines 
a much richer class of models allowing infinite set of commuting conserved operators along with 
their exact eigenvalue solution through the algebraic Bethe ansatz (ABA) |16|| , has not been 
explored. Moreover, the recent controversy around the solvability [|15| as well as the claim and 



denial of the equivalence between some models [H], |TB[ could not be solved satisfactorily due 
to the lack of understanding of the general structures for such systems. 

Recall that the quantum integrability for a lattice model of N sites ensures that, there must 
exist independent and mutually commuting set of conserved operators {Cj},j = 1, . . . , N, the 
Hamiltonian being just one of them. For achieving this one usually shows that a one-parameter 
family of transfer-matrix r(A) = tr(Y\j Lj(X), constructed from the Lax operator Lj(X) and 
generating the set of conserved operators, themselves commute. This in turn follows from 
the Yang-Baxter equation (YBE) satisfied by the associated L,,(A)-operator and the i?-matrix 
solutions. Note that, while in CBA solvable models one can solve the eigenvalue problem 
(EVP) only for the Hamiltonian, the quantum integrable systems allow exact solution of the 
EVP, simultaneously for the whole set of conserved operators, including the Hamiltonian. This 
is done by solving the EVP of the transfer matrix: r(A)|n >= A n (A)|n > through ABA. 
Consequently, for describing quantum integrable systems |L6| one has to start not from the 
Hamiltonian of the model but from the related solutions of the YBE. 

Our aim here is to look into the class of solvable electron and spin models mentioned above, 
from the view point of quantum integrability. We therefore start by discovering the associated 
i?-matrices and the quantum Lax operators for the integrable class in such systems. In doing 
so we exploit certain symmetries of the YBE under twisting and gauge transformations. This 
establishes in one hand the complete quantum integrability of the system and allows to apply 
the ABA method for exactly solving the EVP for the transfer matrix and consequently for all 
the conserved operators, simultaneously. At the same time our formulation helps to identify the 
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integrable subclass of CBA solvable models classified in |15| and find the precise relationship 



between various existing models flTI , [TJ|, 13, [15|] at their Lax operator as well as Hamiltonian 
level. As an important consequence, we detect some error in the derivation of the Zvyagin et 



al model pl| , [12], derive the correct one and thus completely resolve the controversy regarding 
the solvability of the Zvyagin et al model and its equivalence with the Schulz-Shastry model 
[ |l7l , |l8|j . In addition our approach based on the construction of Lax-operators allows us to 
generate new quantum integrable models including generalized t — J models 

The organization of our paper is as follows. Sect. 2. describes briefly the main ideas of 
quantum integrability and ABA. Sect. 3 introduces the symmetries of YBE and the construction 
of the .R-matrix and the Lax-operator. Sect. 4 is devoted to the derivation of models showing 
their mutual relations and establishes the equivalence between different models. The extension 
for our formulation to the multi-chain models and their equivalence are also considered here. 
Sect. 5 generates a new class of models including one for the generalised t — J models. Sect. 
6 presents the exact ABA solutions for the models constructed here. Sect. 7 is the concluding 
section. 



2 Quantum integrability, conserved operators and ABA 

By quantum integrability we will mean the integrability in the Liouville sense. For classical 
models such integrability means the existence of action-angle variables, when conserved quan- 
tities including the Hamiltonian can be expressed through action variables only. For quantum 
models, Liouville integrability similarly demands existence of sufficient number of mutually 
commuting conserved operators {Cj}, the Hamiltonian being one of them. The aim of the 
quantum inverse scattering is to solve the eigenvalue problem for all these conserved quantities, 
simultaneously. Recall in this context that the CBA is designed to solve this problem only for 
the Hamiltonian of the model. Therefore quantum integrable systems are much richer, though 
more involved and demanding than the CBA solvable models and as a rule the quantum integra- 
bility does not follow from its solvability by CBA. For describing integrable systems, naturally 
one can not start from the Hamiltonian as in CBA, but has to introduce some abstract objects 
like the Lax operator Lj(X) and the i?-matrix, which depend on an extra parameter A called 
the spectral parameter . Though the Lax operator is also a matrix similar to the f?-matrix, 
it is a p x p-matrix, while the matrix dimension of R is p 2 x p 2 . Another important difference 
is that, the elements of matrix Lj(X) are quantum operators defined at site j and acting in the 
corresponding Hilbert space, while the elements of the f?-matrix are the usual c-numbers. The 
matrices are considered here in the fundamental representation of the underlying algebra with 
its rank generally coinciding with p — 1. For example for su(2), as in case of the spin-| chains, 
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the Lax operator would be a 2 x 2 and the f?-matrix a 4 x 4 matrix. 

A sufficient condition for the quantum integrability of the system may be given by the 
celebrated YBE 



R ab {X,fi)L a:j {X)L bj {fi) = L bj (fx)L aj (X)R ab (X,fj,), j = 1, 2, . . . , N, (2.1) 
where a, b indicate the matrix or the auxiliary spaces, while j the quantum spaces. The asso- 



ciativity property of the algebraic relation (|2.1| ) yields in turn a similar relation 



R ab (X,^)R ac (X,j)R bc (fi,j) = R bc (fi,~f)R ac (X,~f)R ab (X,fi), (2.2) 

representing a compatibility condition for the i?-matrix. The R(X, /u)-matrix usually depends 
on the difference between spectral parameters A and //, though there may be exceptions as in 
the case of Hubbard model. The well known and well studied integrable models belong to the so 
called ultralocal class, which along with the YBE fl2.1|) obeys also an additional restriction on the 
L-operators: [L ai (X), L b j(fi)} = 0, i ^ j, a ^ b. Since the notion of integrability is intrinsically 
related to the conserved quantities, which are global objects concerning the whole system, for 
its description we must switch over from the local YBE ( |2.1| ) to some global relations. For 
this purpose therefore, we define a global object known as the monodromy matrix, by matrix- 
multiplying the Lax operators: T a (X) = Ylf = i L a j(X), which yields again apxp matrix with 
operator elements T a p. Multiplying the local relations (|2.1| ) successively for all sites j we arrive 
at the global YBE exhibiting the same form: 

R a b(X - fi) T a (X) T b (fi) = T b (fi)T a (X) R ab (X-fi). (2.3) 

Note that the construction of this global YBE is possible only due to the ultralocality property, 
when the Lax operators at different sites can be treated almost like commuting classical objects 
and can be dragged through one another and multiplied to yield ( |2.3|) . This algebraic invariance 
of the tensor product reflects the deep Hopf algebra structure underlying all integrable systems. 
Taking now the trace from both sides of relation (|2.3|) , the i?-matrices get canceled, since 
they can be rotated cyclically under the trace and one gets finally for the transfer matrix 
r(A) = tr(T(X)) the crucial relation [r(A),r(//)] = 0. Defining the conserved quantities as the 
expansion coefficients of 

lnr(X) = J2C n X n , (2.4) 

n 

the commutativity of the set {C n } for different values of n: [C n , C m ] = 0, follows immediately 
and that establishes the complete quantum integrability of the system. 
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2.1 Hamiltonian construction 



Though the starting point in a Quantum integrable system is an abstract Lax operator, it can 
systematically generate the explicit forms of all conserved operators of physical interest includ- 
ing the Hamiltonian of the system. This follows from the fact that the Lax operators through 
successive multiplication leads to the monodromy matrix and the trace of it gives the transfer 
matrix, which in turn using (|2.4|) generates the conserved quantities C n = ^-^lnr{\) \x=o- For 
example, if the Hamiltonian of the model is defined as H = C\ = t'(0)t~ 1 (0) and the Lax 
operator satisfies an important criterion called regularity condition: L a j(0) = P a j, linking it 
with the permutation operator P a j, then the Hamiltonian may be expressed directly through 
the Lax operator as 

H = E ffii+i = 4+i(°) L 7/+i(°) = (2- 5 ) 

3 

Note that due to the regularity condition of the Lax operator one can use the properties of 
the permutation operator: P a jL a j +k = Ljj +k P a j and = tr a P a j = I to derive Hamiltonian 
(|2.5|) , exhibiting nearest neighbour interactions. Fortunately, fundamental integrable models 
in condensed matter physics, e.g. spin models, Hubbard model, t — J model etc. fall into 
the category of ultralocal as well as regular model and their Hamiltonians are with nearest 
neighbour interactions described by C\. The theory of nonultralocal integrable models is still 
in the developmental stage il9j and will not be considered here. 



2.2 Algebraic Bethe ansatz 

As we have stated above, while the CBA is concerned with the exact solution of the EVP for 
the Hamiltonian H, the ABA aims to do so for all conserved operators simultaneously. This is 
achieved by solving the EVP for the transfer matrix, expressed through the diagonal elements 
of the monodromy matrix as r(A) = J2 a T a a{X). The off-diagonal elements of the monodromy 
matrix T Q/3 (Aj) with a < (3, on the other hand correspond usually to the pseudoparticle cre- 
ation operators and their successive actions on the pseudovacuum |0 > generate the n-particle 
eigenstates \n >. The states \n > and consequently the corresponding eigenvalues of r(A) de- 
pends on p — 1 different sets of spectral parameters {A^ },r = 1, . . . ,p— 1 also known as the 
rapidity variables, such that the total number of excitations matches with the pseudoparticle 
number:^ i r = n. Therefore for calculating the related eigenvalue problem, one has to examine 
the actions of the diagonal elements T aa (X)\n >, for which one needs the explicit commuta- 
tion relations between the diagonal and the off-diagonal elements of T(X). This knowledge is 
required for dragging the former type of operators through the later type until they hit the 
pseudovacuum state |0 >. The resulting problem T aa (A)|0 >= a a (X)\0 > can be readily solved, 
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since a a (X) is obtained from the trivial action of the Lax operators on the pseudovacuum state. 
Demanding \n > to be the true eigenstate of r(A), one gets on the other hand some important 
extra conditions, known as the Bethe equations, for determining the rapidity variables {A^ }. 
The exact procedure of ABA depends naturally on the concrete models. Recalling p — 1 as the 
rank of the underlying algebra, we see that for systems with p = 2, e.g. for the spin-| chain 
there is only a single set of rapidities {Aj} , while for m-chain spin models, one gets m sets of 
rapidities {A^ -1 '}, r = 1, . . . , m. Similarly, the Hubbard and the supersymmetric t — J models 
are described by two sets of rapidities {A^° , A^ 1 -*}. Note however that the underlying symmetry 
algebra of the Hubbard model is su(2) x su(2), resulting a 4 x 4 Lax operator, while that for the 
SUSY t— J model is gl(l, 2) with its Lax operator being a 3 x 3 matrix. For models with p > 2 
the ABA steps are recursively repeated giving the nested Bethe ansatz. We shall give some 
details on the twisted Hubbard and the t — J models in sect. 6 for elaborating this process. 



3 Symmetries of YBE and Generating transformations 

For identifying a procedure for generating integrable multi-chain electron models with correlated 
hopping and spin models with interchain interactions we intend to find certain symmetries of 
the YBE, which would yield under particular transformation Lax operators and f?-matrices 
as new solutions of the YBE. Therefore the idea is to start with the L, R matrices of known 
integrable systems and derive the transformed ones satisfying again the YBE and therefore 
representing new integrable systems. From these transformed Lax operators one can derive 
now the Hamiltonians of the desired integrable models with correlated hopping. 
We look for the set of transformations given by 

Rab{X,fi) — > Rab(\,fi) = F ab R ab (X, n)F^\ L a j(X) = F a jL a j(\)Fj a , (3.1) 



and demand that the YBE ( |2.1| ) should remain valid under such a transformation. In what 
follows we will frequently denote L a j by R a j for convenience. It can be checked easily through 



simple algebra that, (34) maps a solution of the YBE into another one, only if the following 
set of relations holds. 

-Rab-F a j Ffoj F b jF a jR abl R a jF b jF ba F^F^jR^^ RfojFfoaFja FfoaFjaRbj^ (3.2) 



along with the conditions that F a ^ must itself be a solution of the YBE (|2.2| ) : 



FabF ac Fbc — F bc F ac F ab . (3.3) 

In our constructions we will explicitly find the transforming operator F by assuming certain 
symmetries. For example, supposing the symmetric condition F ba = F ab = S ab , transformation 
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(|3.1|) looks like a gauge transformation 

R a b(^, y) = S ab R ab (\, n)S~b] L a j(X) = S a jL a j(\)S~^ , (3.4) 
while with anti-symmetry F^ 1 = F ab = T ab , which is known as the twisting transformations 



2(| we get 

R a b(X, fi) = T ab R ab (X, n)T ab ; L a j(X) = T a jL a j(X)T a j. (3.5) 

For model construction we consider transformations with such symmetries only and satisfy 
trivially the essential condition ( |3.3|) by choosing the transforming operators F as mutually 
commuting. Notice that under these conditions all relations in ( |3.2|) can be reduced for both 
the gauge and the twisting transformations to a single convenient form 

RabFacFbc = F bc F ac R atn (3.6) 

where the indices a, b, j are treated in equal footing with R a j to be understood as L a j. 

There exists another type of gauge transformation, where the symmetric operators are given 
in the factorised form: S ab = G ab = g a g b . It can be shown easily that, in such cases no extra 
condition is imposed for the validity of the transformed solutions and the operators g a can be 
arbitrary invertible matrix. To demonstrate this we may focus on the lhs of the YBE (|2.1| ) and 
insert the transformed solutions ( |3.1| ) with the transforming operators in the factorised form to 
get 

lhs = g a gbRab(K fi)ga 1 9b 1 9agjL a j(X)ga 1 9j 1 9bgjL bj (fi)g b 1 gJ 1 

= 9a9b9j(Rab(\y)L aj (X)L b:j (fi))g- 1 g b : 1 gT 1 . (3.7) 

Similar calculations are repeated for the rhs and since the old R, L solutions satisfy the YBE, 
it holds naturally for the transformed solutions. Note that in dragging the factors involving 
c/ a 's out of the original lhs of YBE as done in (|3.7|) , we have used only the mutual cancellation 
of these factors as well as trivial commutativity of the operators acting on different spaces. For 
such transformations as evident from ( |3.7| ) the validity of the YBE condition fl3.3| ) is also not 
required. Therefore the choice of the transforming operators g a can be completely arbitrary. 
We show below that using such freedom we can construct integrable models which can go 
beyond the class of models classified in []I5|. However confining only to mutually commuting 



twisting and gauge transformations, as mentioned above, we can derive an integrable subclass 



of the solvable models considered in [ I5| . We emphasise again that all the models we generate 
here belong to the quantum integrable systems and therefore are much richer than the models, 
which allow only the solution of the Hamiltonian eigenvalue problem through CBA. 
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3.1 Transformed Hamiltonian 



Recall that Q2.5Q links Hamiltonian of a model directly to its Lax operator, if the regular- 
ity condition is fulfilled. We notice further that the regularity condition is preserved under 
transformations (|3~T]), since it gives L aj (0) = F aj L aj (ti)F^ = F aj P aj F~^ = /•'„; /•'„/ / , „ / = P a j- 
Therefore one gets a precise way for obtaining new transformed Hamiltonians explicitly from 



the transformed Lax operators (3.1) as 



3 
3 

= F 33+l ( H 33+l) F jj+l > 



(3. 



where Hjj + \ corresponds to the original model as defined in ( |2.5| ). Therefore if one starts from 
known integrable models like XXZ or Hubbard model, with their Lax operators satisfying 
both ultralocality and regularity conditions, then through (|3.1|) one can derive the transformed 
Lax operator with the same essential properties and in turn can generate through (!T£) the 
Hamiltonian with new interactions introduced by the operator Fjj + \. 



4 Model construction 

For constructing new integrable models following the above scheme we restrict ourselves only 
to the symmetric ( |3.4| ) and antisymmetric ( |3.5|) transformations and build the transforming 
operators out of the Cartan generators H( a ),a = l,...,r, r = p — 1 being the rank of the 
underlying algebra. Due to the mutual commutativity of Hj the transforming operators would 
automatically satisfy the YBE (|3.3|) and therefore the only condition that remains to be satisfied 
for our transformations S a b and T a t> is (|3. . It can be shown that such operators in the general 



case may be constructed in the exponential form |20 



with symmetric or antisymmetric properties on the parameters: 9p a — In Q4.1|) a = 

may also be included in the sum by considering H = 1. Note that the rank r of the associated 
algebra depends on the concrete models and determines also the matrix dimension of the Lax 
operator together with the chosen representation. However since we are interested here in a 
particular class of models with relevance in condensed matter physics, e.g. spin chains, Hubbard 
model, t-J model etc. we would prove the validity of ( |4.1| ) in such particular cases only. 
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4.1 Correlated electron and spin ladder models 

We consider first the Hubbard like correlated electron models involving operators cL^, Cj^,rij^ 
with spin components a = ± and XX Z or XXX type two-chain or the so called ladder models 
involving two independent spin-| operators a, f. Subsequently, we will look into their multi- 
chain generalisations as well as transformed t — J and some other new spin models. Note that 
since through Jordan- Wigner transformation the spin operators can be mapped into fermions 
and vice versa, the Lax operators in case of spin models can be obtained analogous to those 
for the fermionic models, which we consider here in details. Therefore in conformity with ( |4. 1|) 
we propose the symmetric gauge transformations as 

S a j(s) = e ^K(-) n j(+)+ n a(+)™ j( -)). (4.2) 
while the antisymmetric twisting transformation in the form 

7^(0, 7(±)) = e i [ e K(-) n j(+)- n »(+) n j(-))+'>(+)(' 1 a (+)- n j(+))+'T'(-)(' 1 a(-)-^(-)) ? (4.3) 



As we have stated above, we have to check only the validity of ( |3.6| ) for both these operators to 
confirm their usefulness for generating new integrable systems. In all our calculations we need 
only to use the following basic commutation relations for fermions or spin operators like 

f{n a (i3))c] {a) = 5 aj 5 a pJ j(a) f(n a{l3) + 1), f(n a{/3) )c j{a) = 5 aj S af3 c j{a) f(n a{l3) - 1) (4.4) 

[a ± ,a 3 ] = T2a ± , [a VI = ^ (a*) 2 = 0, (4.5) 

etc. We note first that the f? a ft-matrix or the Lax operator of the Hubbard model (simi- 
larly those for the XXZ or XXX spin models) contains strings of operators in the quadratic 
form c a (±\Cb(±)i which are the only parts not commuting with the operators n (±) forming the 
above transformations. However it is evident that, any function of the total number operator 
n a(±) + n b(±) commutes with this quadratic form and hence with the R a b, while the operator 
n c(±) commutes trivially with it. We note now that the transforming matrices in ( |3.6| ) induce 
operators exactly in such combinations: F ac F bc ~ e ico ^*-(K(-)+™K-))" c (+) ± K( + )+n 6(+ ))n c (-)) ) for 

( fOl) and ( |4.3| ), respectively, which based on our above arguments proves the validity of ( |3.6| ) 



and therefore the validity of ( |3.1[ ) for both ( |4.2| ) and Q4.2| ). We may consider also an additional 



transformation in the factorised form choosing 

G aj (g) = g a g 3 , g a (g) = e *f ("-(-)»-(+)). (4.6) 



Since ( |4.2|) , (|4.3| ) and (|4.6j ) can be applied independently, we may consider their combined effect 



by their successive application: 

F a j(0,^(±),s,g) = S aj {s)T aj (9,~f {±) )G aj {g), F aj (0,7 (±) , s, g) = ^/(s)T ai (0, 7 (±) )(G ai ) _1 (^) 

(4.7) 
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for constructing transformation through ( [3.1|) as 

R aj (e,l(±),s,g,\) = F aj (e, li±h s,g)R^ ub \\)F aj (e, 1{±) ,s,g). (4.8) 

Here the transformed R in ( |4.8|) stands for both the changed Lax operator and the i?-matrix 
and naturally satisfies the YBE representing a new quantum integrable system with TV number 
of independent conserved operators. These systems would give a hierarchy of extended Hub- 
bard models with correlated hopping or through Jordan- Wigner mapping from fermion to spin 

operator: n a (±) — > |(1 + ^(i)) anc ^ c i(±) — * °i(±)' c i(±) ~~ ¥ a 7(±)' a hi erarcn y °f s P m ladder 
models with interchain interactions. We have used here the notation ov_\ in place of f for 
convenience. 

The Hamiltonian of such transformed Hubbard systems would be given by ( |3.8| ) with the 



use of the transformation (|4.7| ) as 

H = J2F jj+1 (9 n( ± h s,g)H^F^ +1 (9, 1( ± h s,g) 



j 

= ~ E c ]+i(+) c j(+) /i &!i + c j+n-) c 3H h jjli + Un 3(+) n i(-) + h - c > ( 4 - 9 ) 

3 

exhibiting explicit correlated hopping 

h^l x = e i (- 2 ^( ± )+( ±e - s +3)™ J +i(T)+(± e + s -9)^(=F))_ (4-10) 



Its different parameters keep track of the independent transformations (|4.2|) , (|4.3| ), (|4.6| ) and 
therefore one can study each or any combinations of them by switching off the others. In this 
way, as we see below, one can identify different known models showing also the equivalence 
between them. 



4.2 Connection with general setting and equivalence 



We can now link our model to the correlated hopping electron models classified in [15]. This 
comparison reveals that our model described by the Hamiltonian ( [4.9[ ) with ( 4.10| ) constitutes 



an important integrable subclass of the solvable models of fl5 ], defined by a particular choice 
of their coefficients aj m (cr) and 7j(ct) appearing in the correlated hopping terms as 



ajm(cr) = 5 jm (a9 + s - g) + S j+lm (a9 -(s-g)), jj(a) = -2^. 



where a = ±. Therefore the global unitary transformation defined in |15 

U = exp (iJ2(Zj£ n 3(<*) n H0) + Cj,anj(a))) ■ 



(4.11) 



(4.12) 
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that introduces correlated hopping in the Hamiltonian can also be identified for our case by 
equating the coefficients as 



(T, — fj 
j 



cr,—(T 



1 [s - g), Qj+i = = -o-t 



(4.13) 



with the rest of the coefficients in ( 4.12|) including all being zero. As a consequence all 
coefficients Aj^{a) = 2(£J£ ~ £>j+i,k)i ^ {h 3 + 1) appearing in the models of (l5j also vanish. 
Moreover it is significant to note that the class of quantum integrable models we construct 
correspond to the translational invariant subclass of the solvable models |l5j , as is evident 
from ( |4.11|) , since we have here <x, m (cr) = <x,_ m (cr) with its only nontrivial coefficients «o(±) = 
±9 + s — g and a_i(=F) = ±0 — (s — g), satisfying the necessary condition ao(±) + a_i(=F) = 0. 
Weather this linkage between the quantum integrability and the translational invariance has any 
deeper meaning is yet to be explored. On the other hand we recall that the models introduced 
by Schulz-Shastry Jy|] also exhibits translational invariance, since their coefficients correspond 
to ctj m (±) = =|=<x,_ m with A, ifc (er) = jjf±) = 0. Comparing with (|4.11|) we may notice that, the 
Schulz-Shastry model, in spite of its translational invariance, does not in general correspond 
to our quantum integrable case. However, we will see below that a further restricted model, 
known as the minimal Schulz-Shastry model, given by the particular choice «o — ct-i 7^ 
with all other = does agree with a particular case of the present construction, showing 
that the minimal Schulz-Shastry model belongs to a quantum integrable class. Returning 
again to the present case, we observe that just by restricting to the choice ( |4.11| ) and ( |4.13| ), 
we can use all the formulas derived in [15| for our model and thus can calculate the set of 
boundary phases and other relevant objects for the present case. Moreover, we can find that 
the necessary condition for the CBA solvability of such models prescribed in [IB], holds also 
for our model, since for the choice (|4.11|) it is easy to check that the compatibility condition: 



otj m \cr) + a mj+1 



(— a) — a jm+ i((j) — az m j(— a) = with m — j, j + 1 only, holds true. 



Therefore we conclude that, while all Amico et al model ||15j are only CBA solvable, their 
particular subclass defined by flj.ll|) and (|4.13| ) and generated by our construction not only be- 
comes CBA solvable but also represents an exact quantum integrable system. This integrable 
subclass may be given by the Hamiltonian (|4.9| ) with (|4.10| ) together with all its particular 
cases and is associated with explicit Lax operator and /2-matrix solutions ( |4.8| ). Consequently, 
these models allow a hierarchy of mutually commuting conserved operators having increasingly 
further neighbour interactions, e.g. nearest, next nearest, second next nearest neighbours etc., 
as well as more and more nonlinear interacting terms in their successive Hamiltonians. More- 
over, as we have mentioned already, while CBA aims to solve the eigenvalue problem for the 
Hamiltonian of a single model, quantum integrable systems permit exact eigenvalue solution for 
all its iV-number of conserved operators, simultaneously through ABA. We will demonstrate 
briefly this ABA application to our integrable model (|4.9| ),( ^T0|) in sect. 6. 
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We stress again that one can generate various extensions to the Hubbard model from ( |4.9|) 
by adjusting the independent parameters 7(±), 9, s, g in ( |4.10| ), all of which would be quantum 
integrable systems allowing application of ABA, as well as solvable through CBA models. Their 
Lax operators and i?-matrices would naturally be related to each other by gauge or twisting 
transformations giving a direct way of establishing explicit relationship between the models 
with different types of hopping interactions at their Lax operator level. This throws a new light 
on the question of equivalence of the models, which created significant controversy in recent 
years [ |T8| , |T7| , |T5|] . For example, if we restrict only to the twisting transformation by choosing 
9^0, while switching off the other parameters: g = s = 7(±) = 0, we get an integrable model 
with the Lax operator 

(4.14) 



Raj (9, A) = T aj (6,0)R% u b \\)T ai (6,0) 



with T a j(6, 0) = e t ^ na (-) n J(+>~ ?la (+) n J(->) as in fl4~3| ) and the corresponding Hamiltonian will be 
given as ( |4.9| ) with its hopping interaction reduced accordingly from ( |4.10| ) as 



H 



Z^v c, j+i(+) L, j(+) c T ( -i+i(-) ( -j(-) c 



+ Urij( + )nj(-) + h.c. 



(4.15) 



We observe that ( f4.14p coincides with the L-operator found in [TJj], while ( |4 . 1 5| ) recovers the 



Hamiltonian of the Kundu model [14], which also tallies interestingly with the minimal Schulz- 
Shastry model |T3, I8| . Clearly, from ( |4.15| ) one derives «o — a -i — 9 as the only nonvanishing 
coefficients in its correlated hopping term, which defines the minimal Schulz-Shastry model, as 
mentioned above. We therefore conclude that the minimal Schulz-Shastry model H( mm - ss ) ) 
given by the Hamiltonian (|4.15 ) and identical to the model of [jrj[] represents a quantum inte- 
grable system with its Lax operator and the f?-matrix R^ m ~ SS \9, A) given by fl4.14j) . 

Let us consider now a combination of the gauge and the twisting transformation by adjusting 
the corresponding parameters as s = 9, while switching off the others by choosing g = 7(±) = 0. 
The resulting integrable system will have the Lax operator and the i?-matrix (|4.8|) with the 
transforming operators F a j{9) = S a j(9)T a j(9,0) = e 
e 



2i8(n a 



(-) n j(+)> 



F aj (9) = S-/(9)T aj (9,0) 



" 2i0 ( n a(+) n i(-)) j expressed in the form 



R aj (9,X) = e 2^K ( „ ) n j(+) )) jR (H«6)^ A ^_ 2ie( „, a(+) „ j( _ ) ) 



for the extended Hubbard model and similarly 



R aj {9,X) = e 



-2ie{a 3 a rf) 



(4.16) 



(4.17) 



for the XXZ ladder model. It is exciting to recognise that Lax operator ( |4 . 1 6[ ) (or ( [4.17[ )) is 



exactly same as that for the extended Hubbard model R^(9, A) constructed by Zvyagin et al 
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[ |Tll , |l2j for their single chain case. Comparing ( |4.16|) and (|4.14j) we observe therefore that these 
Lax operators are related only by an unitary transformation and consequently we come to the 
important conclusion that the Zvyagin et al model is a quantum integrable model and is gauge 
equivalent to the minimal Schulz-Shastry model at the Lax operator level, i.e. we have 



(9, A) = S a j(9)R, 



(min—SS) ( 
"aj 



(9,X)S^(9), S a ,(9) 



i[8(n c 



(-)™j(+)+ T M+) n 7(-) 



(4.18) 



Moreover since the Hamiltonian can be determined uniquely from the Lax operator as shown 
in ( |2.5| ), the equivalence between the above two models extends also to their Hamiltonian level. 
However, this equivalence remained obscured [|18| possibly due to the fact that in spite of the 



correct form of the Lax operator in [ 11 , |12| , the corresponding Hamiltonian as presented in the 
same papers JTT], [12| involves some error [21]. As a consequence this form of the Hamiltonian 



appears to be even nonsolvable by CBA, as observed in [|1|, whereas being associated with 
solutions of YBE it should certainly be an exact quantum integrable model! Our approach 
resolves this puzzle by deriving the correct form of the Hamiltonian from ( ]3.8|) using the Lax 
operator R^\9,X) of [|ll], |i~2 |, which coincides with ( |4.16|) . Therefore one concludes that the 
correct Zvyagin et al model should be given by ( |4.9| ) through a reduction of the parameters 
9 = 7(±) = 0, s = 9 in the hopping interaction ( ]4.10|) yielding 



H 



E 



,2ie(n H _ )nj+1(+) )) H (Hub) 



-2j6/(n i( _ ) n J+1(+) )) 



E 4+i(+) c i(+) e H ~ } + 4 + K-) c i(-) e_2 j+iw + Un K+mn + h - c - ( 4 - 19 ) 



It follows therefore from the above arguments that, the proper Hamiltonian of the Zvyagin 
et al model, which should be given by ( |4.19| ) is related through unitary transformation to that 
of the minimal Schulz-Shastry model, which coincides with ( |4 . 1 5| ) . That is one should have 
= UiH^ min ~ ss ^U{ 1 , where the global unitary operator U\ can be constructed easily from 
( [4.12|) and ( |4.13|) as U\ = exp (—i9J2j( n j(+) n j(-))) , which establishes the gauge equivalence 
between these two models also at the Hamiltonian level. An alternative choice of s = —9, leads 
to a complementary model with Hamiltonian (|4.9| ) and correlated hopping /i^j+i = e 2l6n i+ li --'> 
and /iL+i = e~ 2ienj (+), which again is gauge equivalent to the minimal Schulz-Shastry model 
and can be taken also as an alternative but correct form of the Zvyagin et al model. Since 
both of the above models are obtained as particular cases of ( ]4.9|) , they must be CBA solvable 
as well as quantum integrable systems. Based on our explicit results derived below, we may 
draw another important conclusion that, contrary to the Hamiltonian of the models the Bethe 
ansatz results are insensitive to the gauge transformations S ( ]4.2| ) and G ( f4.6| ) and depend only 
on the twisting transformation T ( |4.3|) . Therefore, since at the Lax operator level the Zvyagin 
et al model is S-gauge related to the minimal Schulz-Shastry model, which again is identical 
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to the Kundu model, for all these three models the scattering matrix and the Bethe equations 
related to their eigenvalue problem turn out to be identical. This is explicit from the identical 
Bethe ansatz equations presented in respective works [|ll|, [12], [L3|, [14J and possibly this is the 



reason why the final Bethe ansatz results in [jn], |Tj| turn out to be correct in spite of their 
incorrect form of the Hamiltonian. 

4.3 Multi-chain generalization 

In the above models we have taken a = ± for denoting the |, j components of a s = | single 
chain electron or the two independent spin-| operators in a spin ladder model. For generalizing 
these models to M-chains we have to consider the set of fermionic operators {c]^, c j(a), n j(a)} 
or spin operators a^ a ) extending the range of a to a = 1, ... , M. Note that we can incorporate 
easily multi-chain electron models with spin-| by redefining the indices through a break up 
of a = (7, ±),7 = resulting 4f number of electron chains. Since the rest of the 

arguments go parallelly, we will not distinguish spin-less fermions from fermions with spins in 
what follows, for convenience. Moreover since the idea of construction described above for the 
single chain models allows simple generalization to the multi-chain case, we will give only brief 
account of such constructions skipping the details and stressing only the differences. Firstly 
for generalising the above gauge and twisting operators ( }4.2| ), ( f4.3| ) and ( |4.6| ) to the M-chain 



models, we have to use now more general transformations ( t4.1|) to get the transforming operators 

as 

S^ M ^§) = e ! 2o</3 Sa M n °(°)^(« + ' l °(«' l i(«)) ) (4.20) 
Tlf\d,l,l) = e i ^«</3 e ^^ a(a)njW) " naW)nj(a))+ ^a 7(Q)(na(Q)_nj(Q) ^ (4-21) 

and 

G a f\g) = g a M) g^ M) , g a M \g) = e^Ww^, (4.22) 

respectively, considering the parameters s, 9, g to be symmetric under interchange of indices, i.e. 
0a/3 = Opa etc. The Lax operator and the i?-matrix of the new quantum integrable multi-chain 



model can be constructed generalizing ( [4.8D as 



Rlf\0, 7, 9, A) = F^f\e, 7, s, g)Lif )Hub (X)F^\§, 7, s, g), (4.23) 

where the transforming operators ( [4.7[ ) are generalised by using the combinations of multi-chain 
twisting and gauge transformation (|4.20|) - (|4.22|) 



as 



F^\9,j,s,g) = T^{e,^\s)G^\g), F™{9, 7,^) = T<f>(0, l){S^\s)G^\g))-\ 

(4.24) 
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The set of conserved operators, the Hamiltonian being first in this set, can be constructed in 



explicit form through a straightforward generalization of ( |3.8| ) for the M-chain case by using 
the transforming operators from ( }4.24|) . Since the transformations considered here do not 
affect the interacting term in the Hamiltonian involving Cartan generators, we present only its 
transformed hopping or the XY terms given as 



N 



# (M) = E ^)(^ ! y )( ^ )(i g) ( ^ ; ^ )r l 
3 

N,M 

= - C j+l(a) C j(a)hjj+l + h-C, 

with correlated hopping involving all other chains: 

, (a) _ i (-27(a)+^^^ Q (sgn(a-/3)e a ^(n j+1( ^ ) +n j( ^ ) )+(s a ^-g a ^)(n jW) -n j+1( ^ ) )J 



(4.25) 



(4.26) 



Comparing with the classification of we find that, our Hamiltonian ( f4.25|) is obtained for 
the particular choice of their parameters 



A 



a,/3 
jm 



5 jm (sgn(a - (3)9 a/3 + s af3 - g aj3 ) + 5 j+lm (sgn(a - (3)9 aj3 - (s a/3 - g a p)), jj{a) 



-2 7 (a), 
(4.27) 

with (3 a, which may be considered as the M-chain generalization of the relation ( |4. 1 1| ) . 
Consequently the only nontrivial coefficients appearing in the corresponding global unitary 
transformation (|4.12| ) that introduces correlated hopping in the Hamiltonian would be given 

as generalisation of ([4.13|) . The compatibility of the coefficients ( [4.271) with the necessary 
condition for the CBA solvability [TjJ can be easily checked. Therefore we may conclude that a 
subclass of the CBA solvable multi-chain models of |15| defined through the parameter choice 
( f4.27|) and ([4.28|) should also belong to a much richer class of quantum integrable systems. 
Restricting to the twisting transformation alone, i.e. considering only the nontrivial parameter 
choice 9 ^ 0, we get a model with Hamiltonian ( [4.25D and the correlated hopping term reduced 



from ( [4.26| ) to h^ +1 = e 9 *■ ^ af3< ^ J+1 (' 3 ' + ^W))) _ case j s evidently identical to a 

subclass of the multi-chain Schulz-Shastry model [ 18 1 , known as the multi-chain minimal Schulz- 
Shastry model, defined for the only nontrivial choice of their generating function: A a p(0) = 
A*/3(~l) — 9 aj 3. At the same time, similar to the single chain case, this minimal class tallies 



also with the multi-chain generalization of the Kundu model introduced in ||22|| , where however 
the interchain couplings 9 were restricted only to nearest neighbouring chains. Therefore we see 
that among the multi-chain Schulz-Shastry models only its minimal subclass exhibits quantum 
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integrability and allows sufficient number of mutually commuting higher conserved operators 
together with the Lax operator and the i?-matrix given by ( |4.23| ) transformed by twisting 
transformation ( |4.21| ) with parameters 7 = 0. 



On the other hand considering a combination of twisting (|4.21|) and gauge transformation 



( f4.20p by choosing s = 9, g = one gets the associated Lax operator from ( |4.23p as 



Lif )z (e,x) = C ) (^)^f ) (^o)C )XY (A)(C ) )- 1 W^f ) (^o) 

/ M \ 

= e 2i J2 a </3 e "^ n ^)n j( p) rQL^ )xy (A) ] e ~ 2i ^ a <f )eal3na{l3)nj{a) \ (4.29) 
We notice immediately that the Lax operator (|4.29 ) is a generalisation of ( |4.16 ) and coincides 



exactly with that of the multi-chain Zvyagin et al model |Tl], jl^ . Therefore the Hamiltonian 
generated from it should have the form 

N,M 

ff{M)z = _j- c t +i(a) c. (a)e M(-E /J >„»a^+ic«+Ep<„^w) + h.c. (4.30) 

This also establishes, as in the single chain case, that the correct form of the multi-chain Zvyagin 
et al model should be given by the Hamiltonian ([4.30|) , which is related to the minimal multi- 



chain Schulz-Shastry model by an unitary transformation Um = exp (—i J2j, a ^j3 Oa,pnj(ci) n j(p) j , 
while the Lax operators of these models are connected through the gauge transformation ( f4.20|) 



with s = 6. This therefore resolves completely the controversies raised in all earlier works 
[ |Hj , r8| , [T5| ] concerning the multi-chain case. Restricting to M = 1 one naturally recovers all 



results related to the single chain case obtained above. 

5 Application of the scheme for generating new inte- 
grable models 

Using the symmetries of the YBE, as we have shown above, one can generate new solutions of 
YBE and hence generate new quantum integrable models starting from the old ones. Through 
this scheme we have already identified the integrable subclass of the CBA solvable models 
classified earlier. Using other options of the same scheme we intend to generate now new 
models. Firstly we exploit the freedom in choosing the transformation operator g a ( |4.6j ) in any 
arbitrary form. Note that until now we have chosen these operators from the Cartan generators 
only ignoring their other choices like 

gl = e ld+ ^,g- = e ie ~°7^ (5.1) 

which we may use now in transformation ( |4.6| ) either individually or in combinations. Without 
detailed analysis of individual cases we demonstrate here only a sample case of g + on a nonin- 
teracting two-chain XXX model: H 0xxx = J2j HjJ^ 2 ^ = — J2j + TjTj+i. Using only the 
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CR (|4.5|) it can be easily shown that this transformation introduces various interacting terms 
including interchain hopping as well as pair production terms in the form: 



ttxxx{2) _ V^„+„+ /n-0™(2)w +„- 



Jj+u 

= ^ m - E^O^" W + i - o* +1 ) + a^tf - a?)) - (a^irf - o> +1 ) 
j 

+ ^(ii- ff ?)) + 2 W 7 7feil ( 5 - 2 ) 

We can apply similar transformation involving g~ on top of the Hamiltonian ( |5.2|) to restore 
its hermiticity and generate additional interacting term. Another simpler choice of such gauge 
transformation, given by gf(cr) = e lp±a ^ and similarly for r-operators, act like some rotation 
on the spin operators. Therefore the effect of such transformations on each anisotropic XX Z 
spin chain yields 

Rxxzm = H o xxz{ 2) _ (1 _ A)(*p>+4 +1 + o**+ + x ) + 2(p + ) 1 + (* _> f)). (5.3) 



It is clear that at A = 1, i.e. for the XXX chain the additional interactions in ( |5.3|) vanish 
due the the rotational invariance of the model. Note that in spite of the additional interactions 
appearing in the models like ( |5.2p , ( |5.3| ) they remain quantum integrable by construction and 
the associated Lax operators and R- matrices can be constructed through (|4.7|) using (|4.6|) . Note 
that though such transformations of the spin models were known in one or in the other forms, 
their links with the quantum integrable systems and the solutions of the YBE were perhaps 
never detected earlier. However such gauge transformations, as we have mentioned, do not 
affect the eigenvalue solution and therefore the Bethe equations are given in the same form as 
in the original model. 

5.1 New suppersymmetric t-J models with transformed interac- 
tions 



We can use the freedom of generating operators (fO ) to construct transformed fermionic, 



bosonic or spin models using the higher rank groups. Here we shall consider such an application 
using gl(l,2) for constructing integrable variants of suppersymmetric t— J model with correlated 
hopping and other interactions along with their explicit Lax operator solutions. As we know 
the Id supersymmetric t — J model with the constraint on the coupling constants: J = 2t turns 
out to be an exact quantum integrable system with the i?-matrix and the Lax operator given 
as 

R aj (X) = c(X)I + b(X)P aj , L^\\) = b(X)I + c(X)P aj , (5.4) 

where b(X) = c(A) = ^ and P a j is the graded permutation operator formed by the 
generators of gl(l,2) @. 
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The operators (|5.4j) satisfy now a graded version of the YBE and the regularity condition 



and therefore following the arguments of (|2.5|) the Hamiltonian of this SUSY t — J model 



takes the form H = J2j Pjj+i ■ Since we are interested in the integrable extension of this t-J 
model with correlated hopping, we can apply ( |3 . 1| ) on its Lax operator for generating new 
set of Lax operator and R-matrix solutions of the YBE. We may construct the transforming 
operators S, T, G in the same form (|4.2| ), (|4.3| ) and ( }4.6|) by using pair of operators from the set 
n + , 72,_, uq. However due to the constraint n + + n_ + uq = I we are left with the only choice of 
n± yielding the transformed t — J model as 

h = E(-^EsV) c ^) e_247(CT) + cc -) p 

+ 2t (\s- jS - J+1 k jj+1 + cc + S^ J+1 - in,n i+1 ) + n 3 + n j+1 (5.5) 

where V projects out double occupancies, Sj is the spin operator and rij = nj( + \+nj^ the total 
number of electrons on site j. Note that in ( |5.5|) the operator parts in the correlated hopping 



given by the same form as ( 4.10 ) involving interactions between j and j component electrons 
vanish due to the constraint on double-occupancy indicated by the projector V. However the 
operator parts arise as additional interactions in the spin terms given by 

£ = e -*(2(7(+)-7(-))+E CT =±( e ( rl j( CT )+ n J+i( CT ))+( s -9) CT K'( CT )- n J+i( CT ))). (5.g) 



On top of this transformation another one may be considered in the factorised form ( |4.6| ) taking 
the factors as gf = e lp s j , which is equivalent to a SU{2) rotation. Using the commutation 
relations like 

[S~,Q-\ = [S~,Ql] = 0, [S-, Ql] = Ql [S-,Q+] = -Q-, 

where Q± = (1 — ^•(T)) c i(±)' one can snow that under such rotations, for example with gj, the 
twisted hopping part of the SUSY t — J model fl5.5| ) would be transformed to 

" tV (j2j +1{a) c j{ *)e- 2 ^ + c.c. + ip-(c} +1{+) c jH )(e-^- - e" 2 ^) j V (5.7) 

introducing hopping between up and down electrons. Since [S~,ri( + ) + «(-)] = 0, [S~,n^ — 
Ti(-)] 7^ 0, such transformations preserve the total charge but not the spin, though they have no 
effect on the eigenvalue problem as we have mentioned earlier. Nevertheless such t — J models 
as well as ( |5.5|) with additional interactions represent integrable systems associated with Lax 
operator and i?-matrices obtainable from ( |3.1| ) starting with ( [5.41 ). 

We demonstrate below the application of the ABA method for solving the transformed 
models on the explicit examples of the Hubbard type as well as the SUSY t — J models with 
correlated hopping. 
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6 Algebraic Bethe ansatz for correlated hopping models 



In a quantum integrable system the ABA method, which solves the EVP for the entire set 
of conserved operators including the Hamiltonian by solving the EVP for the transfer matrix, 
works also with equal success for the transformed systems with correlated hopping as introduced 
above. However, as we see below, among all the transformations discussed above only the 
twisting transformations can influence the EVP and yield deformed Bethe equations. Since the 
applications of the ABA to the twisted model goes almost parallelly to that for the original 
models, we present here only the main steps of the ABA formulation for the twisted Hubbard 
flPD and the t-J model flO)) . 



6.1 ABA for the extended Hubbard model 

The Lax operator and the i?-matrix for the Id Hubbard model, which may be given by two 
coupled free fermionic i?-matrices of 6- vertex type, were known for quite some time 

L ( a f b) (Xa) = (L^(X a ) ® {L${X a ))exp{h a a\ a al a ) 
R { a h b ub \X a ,X b ) = (L%\X ab )®(L^\x ab )) cosh h ah cos X ab 

+ (L^\\ ab )®(L[l\\ ab ))(aiy_ a )smhh ab cosX abj (6.1) 

where L^ + > and correspond to 6-vertex free fermionic models and involve two independent 
spin- 1 operators <T(+) and <T(_) respectively, representing two spin- components of the electron. 
In ( |6.1|) the notations h ab = h a — h b , with sinh 2h a = |[7 sin 2X a and X ab = X a — X b , X ab = X a + X a 
denoting the dependence on the difference and the sum of the spectral parameters have been 
used. However, in spite of this fact the explicit ABA formulation of the Hubbard model 
was discovered only recently ||23|| . This happened possibly due to somewhat unusual tensorial 
structure of its 16 x 16 R- matrix and the 4x4 Lax operator. Consequently the monodromy 
matrix has the form 

f B{X) B{X) F(X) \ 
T(X)= C(X) A(X) B*(X) , (6.2) 
^ C(A) C*(X) D(X) ) 

where B(X)(B*(X)) and C(A)(C*(A)) are two component vectors representing one-particle cre- 
ation and anihilation operators, while A(A) is a 2 x 2 matrix with v4 12 , A 2 ± are the corresponding 
spin excitation operators. The scalar operators F(X) and C(A) on the other hand correspond 
to two-particle creation and annihilation respectively. The integrability condition here is given 
by a graded version of the YBE (|2.3|). The ABA becomes tricky due to the presence of a 



number of creation/annihilation operators, which should be treated properly for solving the 
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eigenvalue problem of the transfer matrix r(A) = sTr[T(X)}, generated by the supertrace of the 
monodromy matrix ( |6.2| ): 

2 

[ B W~Y^ A™ (A) +-°(A)]|Ai, . . . , X N ; fit, . . .,fi Nl >= A N>Nl (\, {Xj}, {(ip})\h, ■ ■ ■ , Ajv; A*i> • • • > Vn x > ■ 

a=l 

(6.3) 

This in turn through expansion of lnA^A^A) in powers of A, as seen from (|2.4|) , solves the 
eigenvalue problem simultaneously for all conserved operators including the Hamiltonian of the 
system, e.g. the energy of the model should be given by En = A N)N (0)Aj£ N (0). The most 
important point in solving the eigenvalue problem is to identify a 4 x 4 i?-matrix of the XXX 
spin i-chain hidden inside the .R^^-matrix of the Hubbard model, which turns out to be the 
two-particle scattering matrix S of this system constituted out of the matrix elements of R Hub 
(|6.1|) : Rq q = i?n 5 n, i?7 j7 = i2io,iO) -^10,7 = #7,io- We have to adopt here the nested or repeated 
Bethe ansatz approach, choosing the ferromagnetic vaccuum state as the pseudovacuum |0 > 
of the system, which gives the eigenvalue expression as 

N / D \ N / D 

AWMA.HM) = <0|S(A)|0>n ^(A,^) +<0|D(A)|0>n -^(Mi) 

- n(^(A,A,))A^(A,{A,},{M)) (6-4) 
where the factor related to the second step in the nested ABA is given by 

A^^A^A,},^}) = <0«|Ai(A-A,)|0( 1 )>n |^(A,M 

+ <0«|A 22 (A-A,)|0W>n ^(A,M , (6.5) 

fi V#7,10 / 

where |0^^ > is the nested pseudovacuum with respect to the spin excitations. For identifying 
the effect of the gauge and twisting transfornations we note that we will be concerned only on 
the modifications of the diagonal elements of the monodromy matrix T(X) and the elements 
of the scattering matrix S involved in ( |6.3| ) through ( |6.4j ) and fl6.5| ). It is however not difficult 
to check through such explicit expressions that the gauge transformations G, S do not enter 
in any of these elements and therefore do not have any effect on the eigenvalue formulas and 
consequently on the Bethe ansatz equations. This is the reason why the Zvyagin et al and 
the minimal Schulz-Shastry models related through a gauge transformation exhibit same Bethe 
equations as explained above. We find on the other hand that under the twisting transformation 
( |4.3|) the relevant terms that suffer changes are the transfer matrix operators A 2 2(A — Xj)\0 > 1 = 
e -i29N jjN _ Aj)|0W > and the element of the scattering matrix #7,10 (A) = b(X)e~ 2ie , where 
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6(A) = jrjj- For simplicity we shall consider only the parametr 9 in the twisting transformation 
putting 7(±) = 0. Incorporating these changes in ( |6.5|) we get its modified equation under 
twisting as 



2i8N 



n 



N 



AT, 



Y[ b(X - Xj 



n 6 (a - up) 



(6.6) 



K^-x) 

Notice that at A = Xj and /i/3 singularities appear in some of the terms in the expressions for 
the eigenvalues (|6.4| ),( |6~6l) . Therefore demanding the vanishing of residues at this points we 
can arrive at the Bethe equations for determining the parameters {Xj}, {Hp}- The residue of 
A-N,Ni{X) at A = Xj yields 



A' 



while the same at A 



< 0|S(A i )|0 >= e ip * L 
Hp gives 



e i2eN t 



n 

a=l 



1 H a — Xj + U 



Ha ~ Xj 



J.28N -Q ~ \ + y_ _ -Q ~ + ^ 

j = l ^CK Xj a_y Hot HP 



(6.7) 



(6.1 



representing the Bethe equations for the integrable Hubbard model with correlated hopping. 
Note that for U — > |C7 the Bethe equations (|6.7|) and ( |6.7|) coinsides exactly with the CBA 
results corresponding to the Hubbard model with correlated hopping considered in |T2|, |T^, 



Note again that the system of |TT], |TJ[ is gauge euivalent to the minimal model of ||13|| , which 
coincides also with that of We remind that this equivalence holds at the Lax operator 
and the application of the ABA level, which does not require explicit form of the Hamiltonian. 
The equivalence however extends also to their Hamiltonians, provided they are crrectly derived 
from the Lax operators. 



6.2 ABA of twisted t — J model 

The nested ABA treatment for the twisted Hubbard model presented above, may be carried 
out in a similar way for the twisted t — J model. Interestingly analogous to the Hubbard 
model, we also find here that, though the twisting transformation T as well as the gauge 
transformations G, S considered above contribute in generating the correlated hopping terms 
in the Hamiltonian, only the former type of transformation affect the ABA treatment and 
the related Bethe equations. The essential ingredients needed for the ABA application to this 
system, i.e. the transformed Lax operator and the i?-matrix are given as in (|4.14j) with (|5.4|) 
leading to the corresponding monodromy matrix in the form 

/ An (A) A 12 (X) B 1 (X) \ 



T(X) 



A 21 (A) A 22 (X) B 2 (X) 
\ Ci(A) C 2 (X) D(X) ) 



(6.9) 
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Since the ABA steps for this twisted model is close to those for the original supersymmetric 
t — J model, we will follow the formalism of J/J] pointing out only the essential differences 
under twisting. This supersymmetric model involves the bosonic (B) hole state together with 
the fermionic (F) electron states and we consider the BFF case for concreteness, where in 
the first step of nesting |0 > represent the hole state, over which Ci, C%(Bi, B 2 ) act as the 
creation (annihilation) operators producing fermionic spin- j and j pseudo-particles with charge 
excitations. Spin excitations appear in the next step with ^21(^12) becoming the creation 
(annihilation) operators over the vacuum > mimicing the XXX model. The conserved 
quantities are again given by the supertrace 

sTr[T(X)] = r(A) = -(A n (A) + A 22 (X)) + D(X). (6.10) 

The eigenfunctions are constructed as 

Mo 

|{Af }f°, {A«}* >= II^(Af )|0 > FU MMQ (xf), (6.11) 

3 

where ({A^}) = H^ 1 A X2 (Xp |0 >W F£. Therefore for finding out the eigenvalue we 

have to consider carefully the commutation relations between the elemets of r(A) ( |6.10| ) and 
the creation operators forming the eigenvector (|6.11|) using the matrix relations of the graded 



YBE. This yields finally the eigenvalue expression 

,(0)1 rx(lhx _ Uf \\\L TT 1 A (l) 



M 1 M 

Amo.m^A, {Af }, {A^}) = b(X)) L J] - >(0) AliU(A) + II , M ( ) X ^ (6-12) 



b(X-Xf') V&(Af -A) 

where 

(M Mi 1 M fu\— A (0 M M i 

a£ m (a) = - n fo + (A - Af ) n — m — + n m - n — —m- ( 6 - 13 ) 



Note that the effect of twisting transformation in the above eigenvalue expressions can be 
detected in the terms b±(X) = -^■e ±2ld , where the term b + comes from the transfer matrix 
element, while b- is the R- matrix element R 2 $ appearing also in the two-particle scattering 
matrix. The Bethe equations , as discussed above, can be derived from the analiticity condition 
of the egenvalues yielding 



(6(Af)-W^ = e^n^oT^ 

13=1 Xj - A 'p 



A (0 ) _ A (D + j 



M x (l) x (0) Mi \(1) _ I j 

H x(l) x(0) • - 11 x (l) x (l) „•' 
k=l A« Aj, —I (3=1 A a — 
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Note the similarity between the Bethe equations for the twisted t — J and the Hubbard model. 
However they do differ due to the different structure of their i?-matrices and the choice of 
vacua. The CBA analysis || of the twisted t-J model (|5.5| ) is expected to lead also to the same 
Bethe equations as above, which appeared also in fTl . 



7 Concluding Remarks 

We have identified a class transformations including gauge and twisting transformations by 
exploiting the symmetries of the Yang-Baxter equation, which generates multi-chain integrable 
systems of correlated hopping electron or spin models with interachain interactions. Comparing 



our construction with the classification of such CBA solvable models made in 15], we identify an 



important quantum integrable subclass within this CBA solvable models. The solvable models 



of |L5| , which are not covered by our scheme might still show integrability, though possibly they 



would fall under nonultralocal models JT9| not considered here and need seperate investigation. 
On the other hand our construction can go beyond the models classified in (|TjJ and generate 
other integrable extentions of the t — J and spin models as we have demonstrated here. We have 
applied the ABA method to the Hubbard and the SUSY t — J models with correlated hopping, 
which are integrable models constructed through our scheme. The integrable models allow 
mutually commuting higher conserved operators and the exact solutions of their eigenvalue 
problem through the ABA treatment. Our findings shows explicitly that only the twisting 
transformations need to be incorporated in the ABA equations, while the gauge transformations 
have no effect on them, though both of these transformations can considerably change the form 
of the Hamiltonian. Therefore this draws the important conclusion that all models related 
by gauge transformations share the same ABA equations and the eigenvalues, which is in 
fact the case with the Zvyagin et al model and the minimal model of Schulz and Shastry. 
Since our method is based on the transformation of the Lax operator associated with the 



model and not of the Hamiltonian as done in [15], we can get more general information about 



such transformations. Starting from this general scheme we are not only able to find the 
explicit relationship between different models known in the literature at their Lax operator 
and Hamiltonian level, but also could detect some crucial errors in the derivation of a well 
known model, which resolves completely the recent controvercies around the equivalence and 
solvability of some known models. 

Application of the present scheme of constructing twisted and gauge transformed quantum 
integrable systems to other models of physical interest, like spin ladder ||, t — J ladder fllQf , 



Bariev model |24| should be interesting problems to carry out. Recall also that the transport 



property in strongly correlated electron systems is probed through the Drude weight p5] 
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which is usually calculated using the twisted boundary condition caused by the parameter 
7( + ) = ±7(-) = <fi± in the twisting transformation like ( |4.3| ). Therefore the physical relevance 
of the other twisting parameter 9 considered here, as well as its possible effect in extending the 
notion of the Drude weight should be explored |26| . 

It should be mentioned also that various applications and generalizations of the twisting 
transformation [2(J can be found in [27, 28[] and the references therein. The most relevant 
among them in the present context is [28], where twisting transformations were applied to the 
supersymmetric t — J and U models for deriving their Hamiltonians as well as different forms 
of the Bethe ansatz equations modified by such twisting. 

Acknowledgment The author thanks the AvH Foundation of Germany for its support through 
its Followup Programme and the anonymous referee for pointing out some recent references. 
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